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Abstract 

Phylogeographic studies indicate that many marine invertebrates lacking au- 
tonomous dispersal ability are able to achieve trans-oceanic colonization by rafting 
on buoyant macroalgae. However, less is known about the impact of rafting on on- 
going population-genetic connectivity of intertidal species associated with buoyant 
macroalgae. We hypothesize that such species will have higher levels of population- 
genetic connectivity than those exploiting nonbuoyant substrates such as rock. We 
tested this hypothesis by comparing nuclear multilocus population-genetic struc- 
turing in two sister topshell species, which both have a planktonic larval phase but 
are fairly well segregated by their habitat preference of low-tidal bull-kelp hold- 
fasts versus mid-to-low tidal bare rock. We analyzed population samples from four 
sympatric sites spanning 372 km of the east coast of southern New Zealand. The 
sampled region encompasses a 180 km wide habitat discontinuity and is influ- 
enced by a stable, northward coastal current. The level of connectivity was high in 
both species, and neither of them showed significant correlation between genetic 
and geographic distances. However, a significant negative partial correlation be- 
tween genetic distance and habitat discontinuity was found in the rock-associated 
species, and estimates of migrant movement between sites were somewhat differ- 
ent between the two species, with the kelp-associated species more often yield- 
ing higher estimates across the habitat discontinuity, whereas the rock-associated 
species more often exhibited higher estimates between sites interspersed by rock 
habitats. We conclude that for species with substantial means of autonomous dis- 
persal, the most conspicuous consequence of kelp dwelling maybe enhanced long- 
distance dispersal across habitat discontinuities rather than a general increase of 
gene flow. 



Introduction 

Population connectivity or exchange of genetic material be- 
tween geographical populations is a key topic in marine 
ecology, with a wealth of studies having addressed the role 
of planktonic larvae in facilitating gene flow in otherwise 
nondispersive species (reviewed by Kinlan and Gaines 2003; 
Cowen and Sponaugle 2009). Very few studies, however, have 
addressed the evolutionary importance of marine dispersal 
mechanisms that are linked with habitat association dur- 
ing postlarval life stages. Intertidal ecosystems of the world's 
cool and temperate regions are characterized by large, often 



buoyant, perennial macroalgae that support substantial lev- 
els of marine biodiversity. Whereas some of the invertebrate 
species associated with such macroalgal "forests" (kelp beds) 
are highly mobile and exploit many parts of these three- 
dimensional habitats, many taxa are sessile as adults and 
rely on planktonic larvae for dispersal beyond their home 
site. The common ecological association of sedentary ma- 
rine animals with large buoyant macroalgae may have conse- 
quences for their population-genetic connectivity at various 
spatial scales. Empirical evidence from phylogenetic studies 
of macroalgal-associated species suggests that some taxa as- 
sociated with buoyant macroalgae have been able to disperse 
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well beyond their autonomous dispersal range over fairly re- 
cent (Holocene) times (Donald et al. 2005; Nikula etal. 2010). 
These observations imply that epifauna associated with buoy- 
ant macroalgae could have, owing to passive rafting, greater 
potential for dispersal and gene flow than species inhabiting 
rocky intertidal habitats and/or nonbuoyant macroalgae. 

The large brown seaweed orders Fucales and Laminari- 
ales contain numerous buoyant taxa, including Durvillaea 
antarctica (Fig. 1A), Fucus vesiculosus, Macrocystis pyrifera, 
and Ecklonia radiata (Raven et al. 2005). When such buoy- 
ant seaweeds become detached from their rocky substrate 
by wave action, they and their epifauna may get transported 
over long distances to new locations by currents and winds 
(Fraser et al. 20 1 1 ) . Low levels of population-genetic differen- 
tiation within several buoyant macroalgal species have been 
attributed to gene flow mediated by drifting but still repro- 
ductively active adults (Coleman and Brawley 2005; Muhlin 
et al. 2008; Coleman and Kelaher 2009; Fraser et al. 2010). 
This genetic evidence for transportation of whole adult spec- 
imens between macroalgal stands suggests that drifting of 
macroalgae may also have potential to shape population- 
genetic structuring of "hitchhikers" (i.e., associated epifaunal 
species). 

Rafting is often invoked post hoc to explain unexpectedly 
high levels of long-distance population-genetic connectivity 
in marine coastal species (e.g., Bell 2008; McCormick et al. 
2008; Leese etal. 2010). Nevertheless, rigorous tests of the evo- 
lutionary effects of rafting on dispersal and gene flow of ani- 
mal species associated with macroalgae are generally lacking. 
To control for the effects of other dispersal mechanisms, such 



tests would ideally involve contrasting population-genetic 
structures of codistributed species that share life-history 
characteristics relevant for autonomous dispersal ability, but 
are differentially prone to rafting on macroalgae. In this 
study, we explore the effects of macroalgal rafting on genetic 
connectivity of intertidal invertebrates by studying multilo- 
cus population-genetic structuring of two New Zealand co- 
distributed, phylogenetically sister, intertidal mollusc species 
from the genus Diloma (Gastropoda: Trochidae: Monodon- 
tinae). While a recent genetic study compared large-scale 
genetic connectivity in two other Diloma species, one of 
which is specialized on feeding on decaying fragments of bull- 
kelp in high intertidal and one of which inhabits mudflats 
(Donald et al. 2011), the current study is the first to analyze 
a sister species pair in sympatry to quantify the impact of 
habitat association on population-genetic connectivity. 

The topshells Diloma arida (Finlay 1926) and D. durvil- 
laea (Spencer et al. 2009) are morphologically cryptic species, 
whose speciation is suggested to have taken place as a re- 
sult of habitat specialization (Spencer et al. 2009). Diloma 
arida lives throughout coastal New Zealand in shallow rock 
pools and rocky surfaces immediately above the low inter- 
tidal bull-kelp zone and up into the mid-intertidal zone, 
whereas D. durvillaea is typically found inhabiting holdfasts 
of the southern bull-kelp D. antarctica (Cham.) or, occasion- 
ally, crawling on rock surfaces nearby (Spencer et al. 2009; 
Fig. IB). These separate species were unrecognized until de- 
tailed conchological inspections of a large number of indi- 
viduals found from inside bull-kelp holdfasts and subsequent 
phylogenetic analyses of mitochondrial DNA prompted the 
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taxonomic description of D. durvillaea (Spencer et al. 2009). 
Diloma durvillaea apparently has a more restricted distri- 
bution than D. arida: it has so far been found only from 
scattered locations south of Banks Peninsula on the east coast 
of mainland New Zealand and in the subantarctic Auckland 
Islands (Spencer et al. 2009). Genetic studies of beach-cast 
D. antarctica (Collins et al. 2010) and some of its epifaunal 
invertebrates (Fraser et al. 2011) indicate that transportation 
of bull-kelp rafts takes place at a spatial scale that easily en- 
compasses the whole area where these two Diloma species are 
known to coexist. 

Species of the genus Diloma reproduce by broadcasting 
their eggs and sperm freely into the water column (Grange 
1976). While the mode of larval feeding and duration of 
the larval period have not been specifically studied in the 
genus, studies of many other trochid gastropods (Hickman 
1992), especially of the Monodontinae (Underwood 1974), 
suggest that the larval phase is almost certainly lecithotrophic 
and lasts up to 1 week (cf. Donald et al. 2011). Further- 
more, as both species occupy adjacent intertidal habitats 
that are submerged at high tide when broadcast spawning 
takes place, we postulate that the larval dispersal poten- 
tial of these two species is similar. Differences between the 
species that could impact on population connectivity are 
their different distribution patterns and population sizes: the 
rock-associated D. arida is much more abundant and con- 
tinuously distributed along the east coast of New Zealand 
than D. durvillaea. Assuming these differences have mini- 
mal effect, we can potentially ascribe any major population- 
genetic differences between these two taxa to the effects of 
rafting dispersal of postsettlement D. durvillaea on detached 
bull-kelp. 

More specifically, if effective dispersal by rafting takes place 
in D. durvillaea — in addition to larval dispersal that both D. 
durvillaea and D. arida undergo — we expect D. durvillaea to 
exhibit more frequent migration and hence lower genetic dif- 
ferentiation and genetic distances between populations than 
D. arida (Hypothesis 1). Given that bull-kelp beds, the pri- 
mary habitat of D. durvillaea, may undergo intermittent local 
extinction and regeneration cycles (Donald et al. 2011) and 
also have a patchier distribution than the rocky intertidal sur- 
faces where D. arida occurs, frequent long-distance dispersal 
by rafting may also result in highly unpredictable recruit- 
ment patterns of D. durvillaea, and thus to weak correlation 
between geographic distance and genetic differentiation of 
populations (Hypothesis 2). 

Materials and Methods 

Study area and sample collection 

The large-scale oceanography of our study area in the south- 
eastern coast of New Zealand is characterized by a major 
north-east flowing Southland Current that is generated by 




Figure 2. Sampling sites of Diloma arida and D. durvillaea in the east 
coast of New Zealand. Samples of both species were analyzed for ge- 
netic connectivity from the circled sites; only D. arida samples were 
analyzed from the other sites. Location of the Southland Front bounding 
the Southland Current (direction indicated by arrows) was drawn after 
Chiswell (1996). 

the Southland Front, at the interface between subtropical and 
subantarctic surface waters (Fig. 2). A 180-km-long stretch 
of alluvial mixed sand and gravel beach, devoid of any stable 
intertidal rocks that are necessary for the establishment of 
bull-kelp beds, is located in the northern part of the study 
area along the Canterbury Bight (Eikaas and Hemmingsen 
2006) and separates the northern-most collection site (no. 4) 
from those further south. 

We collected samples of rock-dwelling and bull-kelp- 
dwelling Diloma from eight and five sites, respectively (Fig. 2; 
Table 1), at low tide (Fig. 1 A). Specifically, we obtained suffi- 
ciently large samples of Diloma from both habitats from four 
sites (nos. 1-4) that spanned a total 372 km of coastline: these 
samples provided the basis for between-species comparisons. 
We also collected rock-dwelling Diloma from four additional 
sites, to provide insight into geographic structuring of D. 
arida across distances up to 516 km. In spite of searching, 
we did not find holdfast cavity-dwelling Diloma in these four 
localities, but did collect individuals that were grazing on the 
outer surface of holdfasts and fronds. We collected the kelp- 
associated Diloma (presumed D. durvillaea) from holdfast 
cavities of D. antarctica by detaching the holdfasts from their 
rocky substrate using an axe and/or a wrench bar. We col- 
lected rock-associated Diloma (presumed D. arida) by hand 
from rock surfaces and shallow rock pools in the vicinity of 
bull-kelp beds. We immersed the animals in 96% ethanol 
upon collection, with samples from each habitat stored and 
labeled separately. We replaced the ethanol after 1-2 days, 
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Table 1. Sampling sites, their geographic coordinates, and sample sizes of Diloma arida and D. durvillaea at each site from intertidal rock- and 
bull-kelp habitats. Species identity of each individual was determined by genetic cluster assignment. Samples in parentheses were not used in the 
analyses of population connectivity due to their insufficient size. 







Diloma arida 




Diloma durvillaea 


Site number and locality name 


Lat, long coordinates 


Rock 


Kelp 


Rock 


Kelp 


1 . North of Kaka Point, Southwest of Dunedin 


46° 22.890'S; 169° 47.01 7'E 


30 


1 


2 


16 


2. St Clair, Dunedin 


45° 54.858'S; 170° 29.31 7'E 


31 


0 


8 


35 


3. All Day Bay, Southwest of Oamaru 


45° 11.707'S; 170° 53.850'E 


21 


1 


5 


23 


4. Te Oka Bay, Banks Peninsula 


43° 51.205'S; 172° 47.017'E 


23 


0 


4 


31 






Subtotal 107 




Total 124 




5. Ringaringa, Stewart Island 


46° 54.115' S; 168° 08.803' E 


18 


3 1 


(2) 


(1 1 ) 


6. Jacks Bay, South of Balclutha 


46° 29.947' S; 169° 42.618' E 


28 


0 


(2) 


(0) 


7. South of Kaka Point, Southwest of Dunedin 


46° 24.127' S; 169° 47.185' E 


24 


17 1 


(1) 


(2 1 ) 


8. Shag Point, Northeast of Dunedin 


45° 27.930' S; 170° 49.650' E 


24 


0 


(2) 


(6 + 2 1 ) 






Total 221 









Hound and collected not from holdfast cavities like the other kelp-associated individuals, but from outer surface of bull-kelp. 



and stored the samples at 4°C in the laboratory prior to 
analyses. 

Literature suggests that the habitats of these two species are 
well segregated, with rock-associated D. arida having never 
been encountered inside D. antarctica holdfasts (Spencer et al. 
2009). The holdfast-associated D. durvillaea, on the other 
hand, has occasionally been found to occupy rock surfaces 
adjacent to bull-kelp beds (The Museum of New Zealand 
Te Papa Tongarewa voucher specimens, coll. W. R. B. Oliver 
1920; H.G. Spencer, pers. obs.). Although our observations 
during field collections were generally consistent with this 
reported ecological segregation, we did find some instances 
of ecological overlap: specimens with morphological features 
characteristic of D. arida were occasionally found on the out- 
side of holdfasts and on fronds of D. antarctica. In addition, 
despite the subtle qualitative differences in adult shell fea- 
tures, we note that the two species are not known to have 
readily quantifiable, diagnostic morphological characteris- 
tics, and that subadult individuals can be especially difficult to 
distinguish (Spencer et al. 2009; B. A. Marshall, pers. comm.). 
To discover if any of the analyzed specimens came from atyp- 
ical habitat (i.e., D. arida from bull-kelp holdfasts/D. durvil- 
laea from rock surfaces), and to avoid erroneous inclusion 
of D. durvillaea specimens in the D. arida dataset and vice 
versa, we initially performed a genetic assignment analysis on 
the whole two-species dataset, assuming that two noninter- 
breeding clusters are present in the data, and then used the 
cluster membership estimate of each specimen to assign its 
species identity (see Data analyses for details). 

Microsatellite marker development 

The microsatellite markers for this study were developed 
by constructing an enriched genomic library of D. durvil- 



laea for AG, TG, AAC, AAG, ACT, and ATC repeats fol- 
lowing Glenn and Schable (2005). We screened the library 
for microsatellite-positive clones following Perrin and Roy 
(2000). From the total of 697 clones, we amplified and se- 
quenced 39 clearly positive clones with M13-primers using 
the BigDye Cycle Sequencing Kit and an ABI 3730 Genetic 
Analyser (Applied Biosystems, Foster City, CA, USA). We 
found eleven of the sequenced clones to contain microsatel- 
lite repeats that were shorter than 180 base pairs and therefore 
suitable for visualization on 16 cm long, 8-10% acrylamide 
gels. We designed primer pairs for PCR amplification of these 
microsatellite regions with Primer3Plus (Untergasser et al. 
2007). We checked the loci for the presence of allelic polymor- 
phism by first analyzing six specimens of both species from 
St Clair, Dunedin. We optimized PCR amplification condi- 
tions by working on DNA of approximately 50 specimens 
of each species across all sample localities. After genotyping 
more specimens, we found four of the loci to have a very 
high amplification-failure rate and discarded them. The final 
marker set consisted of seven loci; their primer sequences, 
repeat motifs, amplification conditions, diversity statistics, 
and GenBank accession numbers of the respective cloned se- 
quences are presented in Table SI and were submitted to the 
Molecular Ecology Resources Primer Database. 

Microsatellite genotyping 

We extracted DNA from the topshells using either Chelex 
(Walsh et al. 1991) or an alternative method described by 
Zavodna et al. (2008). We found the latter method to yield 
better-quality DNA template for microsatellite amplification, 
and the majority of extractions were subsequently performed 
using this method. We extracted DNA from a piece of tis- 
sue dissected from under the operculum, a body region we 
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expected to be well preserved due to its immediate contact 
with ethanol. Upon extraction, we retained the dried shell 
and ethanol-preserved soft tissue of every study individual; 
they were deposited in collections of The Museum of New 
Zealand Te Papa Tongarewa in Wellington. 

We amplified DNA of the topshells at seven microsatel- 
lite loci in separate 10-/xl volume PCR reactions. PCR mix- 
ture contained 1-4 fi\ of DNA extraction solution, 0.5 
jaM of each primer, 0.75 U Mango Taq DNA polymerase 
(Bioline, London, United Kingdom), 1 x Tag buffer (Bioline), 
0.8 fiU dNTP, 2.0 mM MgCl 2 and 1 /xg BSA (Invitrogen, 
Carlsbad, CA, USA). We resolved the amplified DNA frag- 
ments on 8-10%, vertical acrylamide gels using the Hoefer™ 
SE 600 Chroma equipment (Hoefer Inc., San Francisco, CA, 
USA), in 10-11 mA current. We ran 10-base pair DNA lad- 
der (Invitrogen) on the edge and middle lanes of each gel to 
size the microsatellite fragments. We stained electrophoresed 
gels in SYBR Green (Invitrogen) solution to visualize the 
fragments, photographed them with GelDoc 2000 (Uvitec, 
Cambrigde, United Kingdom), and later scored the alleles 
from digital and printed photographs without information 
about the population or habitat origin of the individual be- 
ing scored. We controlled the consistency of our allele scoring 
by rerunning some previously scored, heterozygous samples 
alongside newly genotyped samples. Amplification and scor- 
ing success varied between individuals and loci; we restricted 
the final sample set and analyses to individuals that were suc- 
cessfully genotyped at four or more loci. Amplification of 
the DNA samples at more than half of the loci and concur- 
rent amplification failure at some loci despite several PCR 
trials suggests that the missing genotypes were largely a result 
of null alleles (i.e., mismatch of the primers with flanking 
regions). Therefore, we adjusted our analyses for presence 
of null alleles whenever possible. The genotype data used 
in the analyses, complete with habitat- and museum speci- 
men voucher information, were deposited in the Dryad data 
repository (http://datadryad.org/). 

Data analyses 

Initially, we analyzed the total 7-locus dataset, consisting of 
413 individuals of the two morphologically cryptic sister 
species, with the individual-based MCMC-clustering analysis 
in the STRUCTURE software (Pritchard et al. 2000) in order 
to verify/determine the species identity of each individual. 
Assuming that our dataset consists of two biological species, 
we determined the membership proportions ( Q Y , Q 2 ) of each 
individual in two hypothetical populations (K = 2) under the 
model of nonadmixture and independent allele frequencies. 
These analyses ignored our prior expectations about species 
identity based on ecology and/or morphology. We started 
the simulations with a burn-in period of 50,000 iterations, 
followed by 200,000 MCMC iterations. We ensured the con- 
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vergence of the MCMC chains by running five independent 
iterations of each simulation that yielded virtually identical 
results. We deemed individuals receiving a >80% member- 
ship estimate (Q- value) for either of the two hypothetical 
populations to qualify for further analyses as representing D. 
arida/D. durvillaea and discarded individuals with Q-values 
between 0.2 and 0.8 from all subsequent analyses. 

We obtained basic descriptive statistics of the datasets, such 
as total allele numbers and allele size ranges per locus, us- 
ing GenAlEx v. 6.4.1 (Peakall and Smouse 2006). We used 
MICRO-CHECKER v. 2.2.3 (Van Oosterhout et al. 2004) to 
test the loci for Hardy- Weinberg equilibrium (HWE) within 
samples, and to assess whether homozygote excesses in the 
HWE-deviant samples were most likely due to technical arti- 
facts (e.g., null alleles or large allele drop-out) or nonrandom 
interbreeding (e.g., Wahlund-effect). Because of evidence of 
null alleles at every locus, we used FreeNA (Chapuis and 
Estoup 2007) to estimate corrected allele frequencies that 
accounted for their presence. We then analyzed these allele 
frequencies with the software DISPAN (Ota 1993) to ob- 
tain estimates of expected heterozygosity (H s ) within sample 
populations, and total heterozygosity (Ht) that took into 
account the estimated frequencies of null alleles in sample 
populations. We estimated H s and Ht across the total four- 
population sample of each species and from all population 
pairs within species. 

We measured the global genetic structuring of the species 
with three types of statistics — D est , F S j, and G"st — that com- 
plement each other regarding different viewpoints to popu- 
lation structuring, that is population differentiation versus 
gene flow or migration between populations (Meirmans and 
Hedrick 2011). We quantified global genetic differentiation 
among the four sample localities in each species by estimating 
mean D est (lost 2008) over loci, using the package DEMEtics 
version 0.8 (Gerlach et al. 2010) in R version 2.11.1 (R Devel- 
opment Core Team 2010). D est is more appropriate measure 
of population differentiation than _F S t with markers that have 
high allelic diversity, as its maximum value is not limited by 
the within-population heterozygosity as is that of F S t (Jost 
2008; Gerlach et al. 2010; Meirmans and Hedrick 2011). We 
generated 1000 random bootstrap data resamplings to obtain 
95% confidence intervals for mean D est estimates. Either alle- 
les or genotypes were resampled, depending on whether the 
contrasted populations were in HWE at a locus or not, respec- 
tively. To enable estimation of global and between-locality 
migrant numbers (N m ) unaffected by within-population het- 
erozygosity as per Meirmans and Hedrick (2011), we calcu- 
lated global and population-pairwise values of the fixation 
indices F$r and G"st (Meirmans and Hedrick 2011). For 
calculation of Fst, we used FreeNA (Chapuis and Estoup 
2007), which provides estimates corrected for the presence of 
null alleles, and uncorrected estimates calculated following 
Weir (1996). We calculated global and population-pairwise 
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G"st estimates according to Meirmans and Hedrick's (2011) 
Equation 4: 

G" ST = [k-(H T - H s )]/[(k-H T - H s )(l - H s )], 

where total heterozygosity (Hf) and expected heterozygosity 
within sample populations (H s ) were average values over the 
seven loci, obtained with DISPAN as explained above; k is 
the number of sampled populations (k = 4 in case of global 
analyses, and k = 2 in case of population-pairwise analy- 
ses). The estimate of the number of migrants per generation 
(N m ) under the equilibrium island model assumption was 
then calculated following Meirmans and Hedrick (201 1) for 
both species globally and for all population pairs within each 
species as N m = (1 - G" S t)/4Fst- 

We used FreeNA also for calculation of population- 
pairwise Cavalli-Sforza and Edwards' (1967) distances (Dq 
distances) corrected for the presence of null alleles. We 
then analyzed the correlations and partial correlations of 
population-pairwise D c distances, geographic distances, and 
major habitat discontinuity (180-km habitat gap in the Can- 
terbury Bight) by conducting Mantel tests in Isolation By Dis- 
tance Web Service (Jensen et al. 2005). Statistical significance 
levels of the correlations and partial correlations (controlled 
for one of the variables) were tested against 30,000 random 
bootstrap replicates. We obtained the geographic distances 
between the sampling localities by measuring the shortest 
marine route between each locality pair in Google Earth. 

We also applied a fourth analysis method that was recently 
developed for quantification of subtle geographic structur- 
ing especially in species with high gene flow: spatial analysis 
of shared alleles (SAShA) that looks at mean distance be- 
tween co-occurrences of an allele in geographically spread 
data (Kelly et al. 2010). We expected to find a larger observed 
mean distance (OM) in D. durvillaea than in D. arida owing 
to the presumed more frequent rafting potential of the for- 
mer. In addition to comparing the mean values of the statistic 
in the two species, permitted by completely overlapping field 
sampling regime across four sampling localities, we tested 
those datasets and also the eight-population D. arida data 
for the presence of significant geographic structuring of al- 
lele occurrences against 10,000 random bootstrap replicates 
that provided null distributions expected under panmictic 
conditions. 

Results 

Species determination with genetic cluster 
assignment analysis 

We found the expected strong match between multilocus 
genotype cluster and habitat (i.e., D. durvillaea in/on bull- 
kelp versus D. arida on rock) for 76% of the genotyped in- 
dividuals (315 out of 413). However, 13% (or 48) of the 
individuals were strongly assigned to the genetic cluster that 



was unexpected based on their collection habitat (Table 1). 
Twenty-six individuals collected from rock surfaces geneti- 
cally resembled (mean Qi = 0.94, SD= 0.066) those collected 
from inside bull-kelp holdfasts that were preliminarily identi- 
fied as D. durvillaea based on shell characteristics. In contrast, 
only two individuals collected from inside bull-kelp holdfasts 
(at sites nos. 1 and 3) were assigned strongly (Qi > 0.92) to 
the "D. arida" genetic cluster, members of which were pre- 
dominantly collected from rock surfaces and identified mor- 
phologically as D. arida. In addition, 20 individuals collected 
from top of bull-kelp holdfasts or fronds at sites nos. 5 and 7 
grouped genetically (mean Q 2 = 0.98, SD = 0.023) with the 
predominantly rock-dwelling cluster and were subsequently 
regarded as D. arida. Five individuals collected from hold- 
fast outer surfaces associated strongly (Q t > 0.98) with the 
"D. durvillaea" cluster. Eighteen of the "D. durvillaea" clus- 
ter genotypes were not included in further analyses because 
they came from localities where too few "D. durvillaea" had 
been found to warrant meaningful analyses (see Table 1). We 
note that D. arida was found in numbers on bull-kelp only 
in a locality where D. durvillaea was absent from holdfast 
cavities at the time of sampling, and was absent from bull- 
kelp surfaces in localities where D. durvillaea was found in 
holdfast cavities. In summary, our genetic data suggest that 
when sympatric, the frequency at which D. durvillaea occurs 
on rock substrate is far higher than the frequency at which 
D. arida exploits bull-kelp substrate, and that D. arida only 
very rarely utilizes the cavities of bull-kelp holdfasts (Tables 1 
and 4). 

The remaining 12% of the samples — 37 individuals col- 
lected from rock surfaces, three from the blades of bull-kelp, 
and 10 collected from inside bull-kelp holdfasts — were found 
to be genetically "intermediate" (0.20 < Q; < 0.80, SD = 0.19) 
and were omitted from further analyses due to the apparent 
uncertainty of their species identity. Most likely their inter- 
mediate Q-values were due to data missing at the loci that best 
diagnose D. arida and D. durvillaea: 39 of the "intermediate" 
individuals were missing genotypes on at least one of the top 
three between-species differentiated loci (Ddu5-Ddu7) and 
in general had a higher proportion of missing data (23%) 
than the individuals that received Q-values above 0.8 (14%). 
An alternative explanation is that these individuals are inter- 
specific hybrids, although the present dataset is not sufficient 
to resolve this issue. Based on the morphological examination 
of the genetically intermediate specimens by B. A. Marshall 
at The Museum of New Zealand Te Papa Tongarewa, hybrid 
status is not out of the question. 

Polymorphism and null alleles 

Among the total 363 study individuals with clear species assig- 
nation (listed in Table 1 ) , genotypic data were more complete 
in the genetic cluster determined as D. durvillaea (91% of 
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Table 2. Allele size range, number of alleles detected, and range of estimated null allele frequencies at seven microsatellite loci over population 
samples from sites nos. 1-4 of D. arida and D. durvillaea. 



Locus 


Allele size range (bp) 




Number of alleles 




Null allele frequency 




D. arida 


D. durvillaea 


D. arida 


D. durvillaea 


D. arida 


D. durvillaea 


Ddu1 


76-280 


82-280 


15 


17 


0.08-0.23 


0.17-0.22 


Ddu2 


100-270 


112-255 


37 


38 


0.09-0.17 


0.19-0.30 


Ddu3 


155-220 


133-252 


5 


13 


0.09-0.26 


0-0.24 


Ddu4 


108-300 


126-285 


39 


37 


0.19-0.22 


0.10-0.19 


Ddu5 


65-122 


65-140 


9 


16 


0.08-0.23 


0.05-0.09 


Ddu6 


130-152 


133-140 


7 


5 


0-0.13 


0-0.15 


Ddu7 


126-255 


115-197 


41 


21 


0.26-0.30 


0.19-0.21 



single-locus genotypes scored) than in D. arida (83%). Am- 
plification of the locus Ddu7 was particularly unpredictable 
in D. arida, failing in 56% of the individuals that otherwise 
amplified at four or more loci, while its failure rate in D. 
durvillaea was only 16%. Ddu7 was nevertheless included 
in the analyses, because it showed a consistent, albeit low, 
scoring rate across the D. arida localities. 

Allelic polymorphism at the seven microsatellite loci was 
high in both species, with the number of alleles detected 
per locus ranging between five and 41 in D. durvillaea and 
between five and 39 in D. arida (Table 2). Most of the high- 
frequency alleles (frequency > 0.05) were shared between 
species, but high-frequency, species-specific alleles were de- 
tected at all loci except Ddu3. Alleles at the most polymor- 
phic loci (Ddu2, Ddu4, and Ddu7) were typically rare (fre- 
quency < 0.05) and private to one species. 

We observed significant excesses of homozygotes at most 
loci and in most of the sample populations. Nevertheless, loci 
Ddu3, Ddu4, and Ddu5 were usually at HWE in D. durvil- 
laea samples, and loci Ddu5 and DdulO were at HWE in D. 
arida. MICROCHECKER analyses revealed that observations 
of homozygote excess were always compatible with the pat- 
terns expected from presence of null alleles in the data, and 
incompatible with patterns expected from Wahlund-effects 
or allele "drop-out." Null alleles were estimated to be present 
at high frequencies in populations of both species: a summary 
of null allele frequency estimates per sample and per species 
is presented in Table 3. Estimates of null allele frequencies 
per locus in each species are shown in Table 2. 



Table 3. Average estimates and standard deviations (in parentheses) of 
null allele frequencies over seven microsatellite loci in population samples 
from sites nos. 1-4 of D. arida and D. durvillaea. 



Site no. 


D. arida 


D. durvillaea 


1. 


0.171 (0.044) 


0.150 (0.093) 


2. 


0.151 (0.099) 


0.148 (0.097) 


3. 


0.115 (0.097) 


0.192 (0.044) 


4. 


0.158 (0.073) 


0.180 (0.07) 



Subtle genetic structuring but distinct 
geographical patterns of connectivity 

Diloma arida and D. durvillaea both exhibited very low lev- 
els of population structuring overall (Table 4). According to 
D est values and their 95% confidence intervals, both species 
showed low but significant population-genetic differentia- 
tion. Global f S T estimates were not significantly different 
from zero in either species, but this result was expected given 
their high within-population allelic diversity. We did not find 
significant differences in the strength of the overall structur- 
ing between the two species as measured by mean D est and 
F S t, corrected or uncorrected for null alleles, and their 95% 
confidence intervals. The inferred mean number of migrants 
between populations was 1.5 times larger in D. arida than 
in D. durvillaea, but the overlapping confidence intervals of 
the F S t estimates (central for point estimation of N m ) indi- 
cate this difference is nonsignificant. Also, SAShA suggests 
that both species are panmictic over the area of overlapping 
sampling, even though D. durvillaea has approximately 5 km 
wider mean distance between co-occurrence of shared alle- 
les than D. arida. Summed over all population pairs, the D c 
distances were slightly, but not significantly, higher in the 
kelp-associated D. durvillaea (mean: 0.378; SD = 0.032) than 
in the rock-associated D. arida (mean: 0.334; SD = 0.025). 
However, the population-pairwise D c distances of the species 
were strongly negatively correlated; the statistical significance 
level of the negative correlation was 6% (Pearson's correlation 
coefficient r = -0.073, two-tailed P = 0.06; Fig. 3). 

We detected some differing trends between the species in 
the geographic patterns of genetic differentiation. Specifi- 
cally, rock- associated D. arida showed a strongly positive and 
statistically highly significant correlation of geographic dis- 
tances and genetic D c distances when the effect of habitat 
discontinuity was removed (i.e., controlled for; Table 5). In 
contrast, the kelp-associated D. durvillaea showed a nega- 
tive but statistically nonsignificant trend between genetic and 
geographic distances, and that correlation remained non- 
significant even after removal of the effect of habitat discon- 
tinuity (Fig. 4; Table 5). Similarly, the estimates of migrant 
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Table 4. Habitat segregation and global statistics on population-genetic connectivity in D. arida and D. durvillaea, estimated based on microsatellite 
data at seven loci over population samples from sites nos. 1-4. CI — confidence interval from bootstrap-replicating over loci. ENA — estimate obtained 
using null allele frequency estimation. W m — estimated number of migrants exchanged between subpopulations per generation. 





D. durvillaea 


D. arida 


Proportion of the sampled individuals that inhabited bull-kelp holdfast cavities 


85% 


2% 


\A/i+h i n-nrim il atinn hotorn7\/nncih/ W^/.-..., ai/onnorl t~\\/or Inn 
v v I li ill I uuuuid LiUi I r I tr Lcr uz.y y (JjI Ly "SlENA); dvtriaytru Ltvtri iuli 


0 732 


0 737 


Total heterozygosity H T ( ENA ), averaged over loci 


0.749 


0.755 


Mean D est (95% CI) 


0.090 (0.044, 0.135) 


0.083 (0.042, 0.124) 


F ST (95% CI) 


0.004 (-0.003, 0.018) 


0.006 (-0.005, 0.026) 


FST(ENA) (95% CI) 


0.011 (-0.001, 0.033) 


0.007 (-0.002, 0.025) 


G"sT(ENA) 


0.113 


0.118 


W m , as inferred from (I-G^tV^s-rena) 


20.2 


31.5 


Observed mean distance (OM) between co-occurrences of alleles, km 


148.1 


143.4 


Expected mean distance (EM) between co-occurrences of alleles under panmixia, km 


150.0 


146.6 


P-value for the null hypothesis OM = EM, from random permutation test 


P = 0.43 


P = 0.25 



0.37 



. 0.35 4- 
Q 

<u 

c 0.33 1 



o 0.31 -I 
c 

$ 0.29 



r=-0,79 
P = 0.06 



-+ 



-+- 



-+- 



0.32 0.34 0.36 0.38 0.40 0.42 0.44 

Genetic distance: D. durvillaea 

Figure 3. Population-pairwise D c distances of D. durvillaea plotted 
against the spatially coinciding distances of D. arida. Pearson correla- 
tion coefficient (r) and its two-tailed P-value are shown; dashed line 
marks the fitted linear correlation. 



numbers between localities yielded very different patterns 
between the two species. When summed over all pairwise 
comparisons across the Canterbury Bight, the kelp-associated 
D. durvillaea seems to have exchanged more migrants across 



the habitat gap than the rock-associated D. arida (Fig. 5). 
Moreover, we inferred the highest migrant numbers in D. 
durvillaea to have moved between the localities furthest apart 
(between sites nos. 1 and 4) . By contrast, we inferred the rock- 
associated species to have exchanged more migrants than the 
kelp-associated species in two of three comparisons between 
the sites south of the gap, in an area where rocky intertidal 
habitat is never uninterrupted by sandy or gravel beaches 
more than 8 km long, and also between sites nos. 3 and 4, the 
closest site pair separated by the habitat gap (Fig. 5). 

In the extended set of eight D. arida population samples, 
the correlation between genetic and geographic distances was 
positive, but not significant at the conventional alpha-level 
0.05. With the effect of habitat discontinuity removed, the 
correlation was not significant either (Table 5). SAShA of the 
wider D. arida data yielded a 1 38.3 km mean distance between 
co-occurrences of alleles, which was not significantly different 
from the expected mean distance under panmixia, 141.3 km 
(P = 0.28). 

Discussion 

In order to obtain a better understanding of the mechanisms 
and spatial scales of population-genetic connectivity in the 



Table 5. Coefficients of correlation and partial correlation between Cavalli-Sforza and Edwards' genetic distances corrected for null alleles (Gen), 
geographic distance (Geo), and major habitat discontinuity (Gap) between population samples of D. durvillaea and D. arida, and their statistical 
significance level from Mantel test (one-sided P-value). Pairwise comparisons across the Canterbury Bight (see Fig. 2) were classified as involving a 
major habitat discontinuity. The results on the first row are presented graphically in Figure 3. 



Variables 


D. durvillaea sites nos. 1-4 


D. arida sites nos. 1-4 


D. arida sites nos. 1-8 


Gen, Geo 


-0.39 (P = 0.76) 


0.50 (P = 0.25) 


0.27 (P = 0.22) 


Gen, Geo; controlled for Gap 


-0.53 (P = 0.76) 


0.94 (P < 0.00) 


0.27 (P = 0.21) 


Gen, Gap 


0.51 (P = 0.33) 


-0.77 (P = 0.08) 


0.18 (P = 0.08) 


Gen, Gap; controlled for Geo 


0.61 (P = 0.76) 


-0.97 (P < 0.00) 


0.20 (P = 0.07) 
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Figure 4. Pairwise Dc distances of four codistributed D. durvillaea and D. arida populations, and of eight D. arida populations, plotted against 
geographic distance between the sampling sites. Closed circles — site pairs connected by rocky intertidal habitat; open circles — site pairs separated by a 
180-km habitat discontinuity at the Canterbury Bight. Dashed line — reduced major axis regression for the total correlation of genetic and geographic 
distances, uncontrolled for habitat discontinuity. Global and partial correlations and their Mantel test results are reported in Table 5. 



Diloma durvillaea ^^^^^^^ 

y \ 




Figure 5. Estimates of migrant numbers (W m ) exchanged between pop- 
ulations of D. arida and D. durvillaea during a generation, calculated as 
A/ m = (1- G"sttena])/(4Fst[ena]). °°— indefinite, owing to F SJ < 0. 

marine environment, Kinlan and Gaines (2003) called for 
more attention to be paid to interactions between species 
that make up a community. Our study addressed the effects 
of interaction between a habitat- forming species of an inter- 
tidal community (D. antarctica) and grazing gastropods (D. 
durvillaea, D. arida), which differ in their reliance on and 
occupation of D. antarctica beds. 

Long-distance gene flow appears higher 
in macroalgal-associated species 

Contrary to predictions, we did not find differences in the 
overall level of genetic structuring or in the magnitude of 
genetic distances between populations between the rock- 



associated D. arida versus kelp-associated D. durvillaea (Hy- 
pothesis 1). Both species apparently maintain high levels of 
gene flow across the study area. Nevertheless, in keeping 
with the predictions made at the outset of this study, we 
detected somewhat differing geographical patterns of con- 
nectivity (migrant numbers; correlations of geographic and 
genetic distances) between the two species (Hypothesis 2), 
suggesting that habitat associations can have some evolu- 
tionary consequences for intertidal marine taxa at the studied 
spatial scale. Specifically, populations of D. durvillaea were 
overall more connected across the major Canterbury Bight 
habitat discontinuity than were populations of D. arida. By 
necessity, D. arida must rely on larval dispersal alone to tra- 
verse this 180 km wide habitat gap, whereas D. durvillaea has 
the added potential to use rafting as its dispersal mechanism. 

Our results indicate that the Canterbury Bight represents 
a weaker barrier to gene flow of D. durvillaea than of D. 
arida. The high migration rate estimates of D. durvillaea 
across the habitat discontinuity are consistent with recent 
studies of beach-cast bull-kelp in this region. In particular, a 
genetic survey showed that oceanographic transport of de- 
tached bull-kelp drifting across the Canterbury Bight (e.g., 
between study sites nos. 2-4 and 1) occurs on a regular ba- 
sis (Collins et al. 2010),which suggests that higher migration 
rates of D. durvillea across the Canterbury Bight from sites 
nos. 1-2 could indeed be due to rafting events. The observed 
high levels of allelic diversity at site no. 4 at all studied loci 
suggest that gene flow across the Canterbury Bight is ongoing 
(as opposed to, e.g., a single colonization event). Rafting dis- 
persal may be particularly important for local recolonization 
and population recovery in D. durvillaea, considering the 
susceptibility of D. antarctica to local population extinctions 
(Donald etal. 2011). 
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Based on our comparison of genetic connectivity of a kelp- 
associated and a rock-associated intertidal snail, we expect 
that species inhabiting strongly buoyant macroalgae that 
grow in environments affected by strong tidal and wind- 
driven currents are predisposed to gene flow by rafting dis- 
persal. Rafting dispersal appears to have potential to enhance 
connectivity between populations that are not within easy 
reach of dispersal achieved by planktonic larvae. The absence 
of genetic isolation-by-distance/habitat gap pattern in the 
kelp-associated D. durvillaea and its presence in the rock- 
associated D. arida — albeit based on a very small number 
of datapoints and only one major habitat gap — implies that 
gene flow effects of rafting dispersal can be highly unpre- 
dictable; in particular, they are not necessarily dependent on 
geographic distance even under influence of a stable, unidi- 
rectional major current. 

Our study focused on a kelp-associated snail that exploits 
only a very small part of the alga, the holdfast. We believe 
holdfast-associated species are particularly likely to survive 
rafting journeys, as the complex structure of holdfast cavi- 
ties and tunnels of many large seaweed species provides good 
shelter from wave action. In spite of constituting only a small 
proportion of the total biomass of typical large macroal- 
gae, holdfasts support far more epifaunal biodiversity than 
blades or stipes (Christie et al. 2003; Wlodarska-Kowalczuk 
et al. 2009). Therefore, the microevolutionary effects of raft- 
ing dispersal potentially affect a considerable part of marine 
intertidal biodiversity. Evidently, the scale at which rafting 
dispersal impacts population connectivity of a species is de- 
pendent on the species' pelagic larval duration and maximum 
larval dispersal distance. Evolutionary effects of rafting dis- 
persal are likely to be more pronounced in kelp-associated 
species that lack planktonic larvae, unlike D. durvillaea. 



Local gene flow seems higher in the 
rock-associated species 

In line with the habitat preferences of the two snail species, 
our study suggests that D. arida has exchanged (on average) 
higher numbers of migrants across local scales in southern 
South Island (e.g., between sites nos. 1 and 3) where no sandy 
or gravel beaches longer than 8 km separate rocky intertidal 
habitats. As D. arida is abundant across this rocky region 
(with many additional populations that were not sampled 
for this study), there are enhanced opportunities for local 
gene flow. By contrast, local gene flow is presumably weaker 
in D. durvillaea because of the more limited bull-kelp habitat 
whose presence is more strictly dictated by the slope and wave 
exposure of the shore. The very large effective population 
size of D. arida (relative to the rarer and more restricted 
D, durvillaea; Spencer et al. 2009) may also partly account 
for the apparent lack of strong genetic differentiation among 
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sampled populations, due to weaker effects of genetic drift 
(cf. Donald etal. 2011). 

Despite the uniformly directional (south-to-north) flow of 
the Southland Current (Chiswell 2009), there is a possibility 
that variable rates of flow in this region could partly explain 
regional variation in gene flow estimates. For instance, the 
surface speeds recorded between sites nos. 1 and 2 have been 
estimated to be five times slower (typically below 10 cm/s) 
than in the region north of the Otago Peninsula between 
sites nos. 2 and 4 (Chiswell 1996). Thus, the relatively high 
speed of the Southland Current in the Canterbury Bight area 
could perhaps have contributed to the higher estimate of D. 
arida migrants between sites nos. 3 and 4. For example, at a 
steady speed of 50 cm/s, larval offspring from site no. 3 could 
potentially cross the habitat gap and reach site no. 4 within 
5 days — a time interval that is roughly concordant with the 
life span of trochid veligers (Hickman 1992). 

The validity of the inferences we have made on the pat- 
terns and strength of population-genetic structuring in the 
two Diloma species depends on the comparability of the mi- 
crosatellite data between the species. Given that we inferred 
null alleles to be abundant in the data, it is important to con- 
sider whether they could have biased our between-species 
comparisons of patterns and strengths of population-genetic 
structuring, on one hand by disproportionately inflating all 
population-pairwise F S t and D c estimates (Chapuis and 
Estoup 2007) of one species relative to the other, and/or on 
the other hand by influencing different population-pairwise 
contrast in one species than in the other. Given that the esti- 
mated null allele frequencies were fairly similar between the 
species across populations and loci, we contend that the null 
allele correction methods of Chapuis and Estoup (2007) — 
that decrease F S j and D c — were equally efficient for the two 
species and that null alleles have therefore not seriously biased 
our between-species comparisons. We also tested to see if the 
locus Ddu7 (whose amplification success differed greatly be- 
tween the species) had affected the between-species compar- 
isons and found out that it had not: for example, population- 
pairwise Dc distances from analyses that omitted Ddu7 cor- 
related closely with those from the seven-locus data in both 
species (Pearson r > 0.98), and the absolute differences in 
the distances were small (D. durvillaea, mean AD C = 0.016, 
SD = 0.006; D. arida, mean AD C = 0.018, SD = 0.009). 

Our study suggests that for intertidal snails with a plank- 
tonic larval dispersal stage, ecological association with a 
buoyant macroalga does not confer dramatic consequences 
to population-genetic structuring, at least on the geographic 
scale and genetic resolution achieved by our study. However, 
given that D, durvillaea has a very fragmented and sparse 
distribution compared to D. arida, it is remarkable that the 
two species exhibited comparable levels of genetic connectiv- 
ity. Either the less common D. durvillaea populations must 
produce more gametes and larvae than D. arida, or their gene 
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flow is enhanced by some other mechanism such as rafting of 
reproductively active adults. Genetic studies of species pairs 
with more balanced population distributions and population 
sizes are still needed to more robustly test the evolutionary 
significance of rafting dispersal for epifauna associated with 
buoyant macroalgae. To this end, we are preparing a compari- 
son of population-genetic structuring in rock-associated ver- 
sus kelp-associated chitons of the genus Sypharochiton along 
the eastern coast of New Zealand, using AFLP-markers. Im- 
portantly, this additional study will provide much needed 
phylogenetic replication for testing the dispersal-by-rafting 
hypothesis. 
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